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Loop quantum cosmology allows to replace the original singularity by a quantum bounce followed 
by a phase of fast expansion called super-inflation. In this paper, we use a simple analytic description 
of super- inflation using the equation of state of a fluid. We then derive the power spectrum of tensor 
perturbations produced. The effects on the B-mode polarization of the cosmological background 
radiation are discussed. For a large region of phase space, the spectrum should be distinguishable 
from the standard inflationary prediction for future observations. 

PACS numbers: 



I. INTRODUCTION 

The Loop Quantum Cosmology (LQC) [J, [1] (a re- 
cent review) is mainly the application to homogeneous 
spacetimes of the quantization methods of Loop Quan- 
tum Gravity (LQG) [3HS], a non-perturbative and back- 
ground independent canonical quantization of general rel- 
ativity. A central result of the theory is that the singular 
evolution encountered in the Wheeler-de Witt theory no 
longer appears [SI [7] . This is commonly interpreted as 
the replacement of the big bang singularity by a bounce. 
That result is a strong motivation to obtain quantitative 
predictions from LQC as a cosmological scenario. Since a 
complete quantum solution is not available at present, a 
constructive means to obtain some insights of the conse- 
quences of the phase close to the bounce is that of using 
effective equations [H H] . In such equations, the effects 
induced by the quantum nature of the geometry are taken 
into account through corrective terms in classical equa- 
tions. In isotropic flat models (which are the only models 
considered here) two main corrections are used : the " in- 
verse volume correction" which comes from the spectra 
of the quantum operator related to inverse triads, and 
the " holonomy correction" coming from the quantization 
scheme in which holonomy of the connection are basic 
variables. The relative importance of the two type of 
corrections is still under discussion, it is then usual to 
consider their effect separately. In the present work we 
focus on holonomy corrections and more specifically on 
the possible effects of the super-inflationary phase (in 
which the Hubble factor is growing) [TOl E] on the ten- 
sor modes of the power spectrum of the CMB. A number 
of study have already been devoted to tensor mode in 
LQC [m . Recently, effect of the big bounce on the 
B-modes of the CMB has been considered j23j (see also 
[M] for a different approach), where the matter content 



of the universe was taken into account through a massive 
scalar field. In the present work we modclizc the con- 
tent of the universe during the phase of super-inflation 
using the simple equation of state p — ujp, for which the 
modified Friedmann equation admits solutions in closed 
form for some relevant values of the constant uj. We show 
that during the super-inflation that equation of state for 
a massless scalar field (w = 1) is a very good approxi- 
mation to the massive scalar field in a fast roll regime 
(0^ >> V{(j))). We compute in both cases (massless 
scalar and fast rolling massive scalar) the power spectrum 
after the phase of super-inflation including the holonomy 
corrections. We then estimate the B-modes and discuss 
the sensitivity of the spectrum to the model considered. 

In Sec. |Tl]we describe the dynamics the of a fluid gov- 
erned by a simple equation of state p = tup. We study 
how this equation compares with that of a massive scalar 
during the super-inflationary phase. The primordial and 
B-mode power spectra for the case uj = 1 are determined 
and discussed and compared to the standard inflationary 
prediction in Sec. 
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Except otherwise stated, the ge- 
ometric units h = c = G = I are used throughout the 
paper. 



II. DYNAMICS OF THE SUPER-INFLATION 
FROM A SIMPLE FLUID MODEL 

We assume that the dynamics of the universe is de- 
scribed by the LQC modified Friedmann equation [3]: 
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where n — SttG and pc is the so-called critical density. 
The modified Friedmann equation ([T]) captures the ul- 
traviolet corrections arising from the discrete quantum 
geometric effects. In particular, we see that when p = pc 
the theory departs drastically from general relativity: the 
Hubble parameter vanishes, and the Universe experiences 
a turnaround in the scale factor. The Big Bang singu- 
larity is therefore avoided and replaced by a big bounce. 
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The critical density can be expressed as a function of the 
so-called Barbero-Immirizi parameter 7 [2]: 
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Combining ([l]) with the derivative of the continuity 
equation p + H[p + p) = leads to 
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which shows that for p ^ the Hubble factor increases 
between the bounce (p ~ pc) and p = Pc/'^- This is the 
phase of super-inflation. 

As for most non linear equations, it is often impossi- 
ble to find an exact solution of the modified Friedmann 
equation. However, with the generic equation of state 
p = up, which, using the equation for the energy conser- 
vation, translates in p = ^a"^, where ^ and X = 3(1 -t-w) 
are constants and a is the scale factor, the modified Fried- 
mann equation becomes 
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This equation can be solved literally for certain values of 
X, in particular for X = (inflation or cosmological con- 
stant), X = 2 (unknown fluid), X = 4 (ultra- relativistic 
matter) and also (not so easily) for X — 6 (massless 
scalar fleld). 

The boundary condition for the modifled Friedmann 
equation is given at the bounce, where the density reaches 
the critical density. Thus ^ = PcOff, where ao is (chosen 
to be) the first eigenvalue of the scale factor which, ac- 
cording to [55], is: 
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Note that, we can also determine thanks to ([s]) the value 
of the scale factor at the end of super- inflation (p — Pc/2) 
for X >0: 

a-end = 2^ao, 

thus, the number of e-folds of super-inflation is simply 
given by: 



N : 



ln(2) 
X 



It is then clear that for X > 1 super-inflation alone won't 
give a decent amount of efolds. 

In order to compare with the the case of a single mas- 
sive scalar field one can rewrite the equation of state for 
the scalar field under the form p = up with uj and thus 



X time dependant. Since the evolution of the scalar field 
is as usual governed by the Klein- Gordon equation 



3H(f) + m^(l) = 0. 
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By restricting the study to the case where the back- 
reaction is negligible |26j . one can solve numerically this 
equation together with the modified Friedmann eq. ([T]) 
to obtain the field (j){t). Then, from the expression of the 
pressure and the density of a scalar field, namely 



one obtains 



X{t) = 



(5) 



which is here a function of the time which depends on 
two physical parameters: m and Fb = m^0(O)^/2pc- It 
is pictured in Fig. [l] 

The evolution of X(t) shows three main phases: a pre- 
inflationary phase {Q < t < 10'^) during which w 6, 
an intermediate phase (10"^ < t < 2.10^) during which 
X(t) decreases and flnally a phase of slow-roll inflation. 
The length of the pre-inflationary phase is related to the 
value of Fb- if the potential energy at the bounce is 
small, it will take longer to reach a slow rolling regime. 
During this pre-inflationary phase the dynamic of the 
massive scalar field can be approximated by the simple 
fiuid p — up with UJ ^ 1 {X ~ 6). Such an equation 
of state corresponds to a massless scalar field. A more 
accurate description of this phase can be obtained with a 
fast rolling scalar field, i.e. such that 0^ >> 1^(0), with 
equation of state X{t) k, Q. In what follows we will study 
these two approximations. Note that the time tend at 
which the super-inflation ends is easily determined from 
p{t) (Fig. 2]) thanks to the relation p(iond) = Pc/2. One 
observes that the super-inflation represents a very small 
fraction of the pre-inflationary phase. 
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FIG. 1: Evolution of X{t) for m = 6.10"^ and Fb = 10" 
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its solution is then: 
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where d is a constant determined by the initial conditions. 
At t = 0, one gets 



Oo = — = — (1 + Kpctf), 
Pc Pc 



(12) 



which implies that d — 0. The conformal time is then 
given by 



FIG. 2: Evolution of p{t) for m = 6.10"^ and Fb = 10"^ The 
horizontal dashed line represents the value at which super- 
inflation ends. Here super-inflation lasts a bit more than 0.1 
and thus 10* times less than the phase with X ~ 6. 



Let us now turn to the tensor modes. From the consid- 
erations of the previous paragraph, one can assume that 
the matter is described by a scalar and then consider the 
two approximations (massless and fast roll) for the pre- 
inflationary phase. In order to obtain the power spectra 
of the perturbations, we have to solve the equation of 
propagation of the tensor modes, which can be written 
as 
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Ukiv) = 0, (6) 



where ra\^i^ is the holonomy correction |25j 

mL.(^)=2Ka2(r;)i^Qp-y(0)). (7) 

Thanks to ^,V{(j)) can be written as function of X: 

X 



V = ^a-^ 1 - 
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and thus 
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We now determine the two quantities rnf^^i^ and 
a" (77) / 0(77) for each approximation: free massless and fast 
rolling scalar. 



dt 
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where 2^i("-; d, z) is a hypergeometric function. 

One can compute the terms of the propagation equa- 
tion: 
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and 
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B. Fast rolling scalar field {X ^ 6) 

For a fast rolling scalar field (^0^ >> V{(j))), the F 
man ([T]) and Klein Gordon Q equations reduce to: 



iH(t) = 0. 



(16) 
(17) 



Thus 



A. Free Massless scalar field (X=6) 

The modified Friedmann equation can be written as: 



whose solution is 



(18) 
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where A is a constant fixed by the initial conditions. At 
t = 0, p = pc = \(tP and so one gets 



<^(0) 




(20) 



which implies that ^ = 0. Using the Friedman equation, 
one can then derive the Hubble parameter: 



H{t) 



1 npct 
6 3 + Kpct"^ 



and since H — a/a, the scale parameter is: 



(21) 



where rjend corresponds to the end of the pre-inflationary 
phase (tend ~ 10^, see Fig. [l] ). Thus, we are not 
only computing the pertubations produced during super- 
inflation but also those emerging from the intermediate 
phase that takes place just before the slow-roll regime. 

In both cases (Fig. [s] and |4]) , the slope at large scales 
grows like and the spectra begin to oscillate very 
slightly around k = 10~^. The minor differences be- 
tween the two spectra can be ignored (see Fig. |5| as they 
would be undistinguishable at the scale of the CMB. We 
therefore conclude that a simple massless scalar field is 
suitable to describe the perturbations produced by super- 
inflation. 
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As in the previous case, the conformal time takes the 
form 
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We can compute the terms of the propagation equation: 
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FIG. 3; Power spectrum of the pre-inflationary phase mod- 
elized as a massless scalar field (X — 6). 



and 



33 



upcO-o ^npcfi — 9 



(3 + upct^) 3 
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which are, contrary to the conformal time, slighty differ- 
ent from those obtained for the massless scalar field. 



III. POWER SPECTRA 

A. super-inflation and beyond 

Let us consider the power spectra of the super-inflation 
alone. For both our modelisations, the equation of prop- 
agation of the tensor modes admits no exact solution 
and was thus solved numerically using a Runge-Kutta 
method, the power spectra could then be expressed as: 



27r2 
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(26) 




FIG. 4: Power spectrum of the pre-inflationary phase mod- 
elized as a fast rolling scalar field [X « 6). 
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B. Complete spectrum of perturbations 

To obtain the power spectrum corresponding to the 
total primordial epoch we have to connect the pertur- 
bations that arise from the pre-inflationary phase to a 
phase of standard slow-roll inflation. To this end we as- 
sume that we neglect the perturbations produced during 
the transition between the two main phases. We also con- 
sider that the quantum corrections are negligible during 
inflation since p « pc [26]. The perturbations produced 
during inflation thus take the usual form 



Vkiv) = -CJ^ Hi'\-krj) -Cj"^ Hi'\-krj), 

y 4 2 y 4 2 

(27) 

where Ci and C2 are constants and hI^' and Hankel 

2 2 

functions of the first and second kind. At the end of the 
first phase, we impose the continuity of the perturbations 
and of their first derivative: 

UkiVend) = Vk{r]end), 

u'kiVend) = v'kiVend)- (28) 

We can then produce the power spectrum of Vk at the 
"end" of inflation (so that the whole phase would yield 
70 e-foldings). As in the case of pre-inflationary phase 
alone, the power spectra (Fig. [6] and [7]) still behave like 
at large scales. For large k the spectra are nearly 
scale invariant as one would expect with inflation. Some 
oscillations appear between those two regimes. Lowering 
the value of the parameter Fg increases the amplitude of 
these oscillations which thus appear as the consequence 
of the bounce. As seen in Sec. decreasing Fb increases 
the length of the pre-inflationary phase, thus giving more 
time to the holonomy correction to modify the spectrum 
of perturbations. 
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FIG. 7; Power spectrum of the pre-inflationary phase followed 
by inflation for m = 6.10~^ and Fb = 10~^^. As usual only 
the super-Hubble modes are considered. 



C. B-mode angular power spectrum 

Since we have computed the power spectrum of the 
tensor perturbations, we are now able to produce the 
primordial part of the B-modc angular power spectrum 
of polarization of the CMB. Wc do not know the modifi- 
cation of the spectrum of scalar perturbations produced 
by the holonomy correction, we will thus assume that 
scalar perturbations are given by the usual nearly scale 
invariant spectrum. This will allow us to produce the 
part of the spectrum that emerges from the conversion 
from E-modes to B-niodes through lensing. Moreover, 
the holonomy correction is an ultraviolet correction while 
the lensed part of the BB angular spectrum is supposed 
to dominate mostly at small angular scales. 

The B-mode power spectra are obtained using CAMB 
|27) . For the cosmological parameters, we took those 
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FIG. 8: Parametriz atio n of the primordial power spectrum 
given by equation pg]) with (fco,i?) = (10"^500) (solid), 
(fco,i?) = (10-^ 500) (dashed) and (fco,i?) = (10-^5000) 
(dotted). 



FIG. 9: Angular power spectrum of the B modes as a function 
of fco with a tensor-scalar ratio of 10~^ and R = 500. 



given for kCDM by WMAP+BAO+HO [H]. As in- 
put for the primordial power spectrum we used a simple 
parametrization of the previous spectrum inspired by |23| 



P(fc) 



R 



(29) 



The parameters /cq and R respectively determine the 
position and amplitude of the transition between two 
regimes (Fig. [8|: P[k) oc at large scale and the spec- 
trum is scale invariant at lower scales. The transition 
between those two regimes presents a bump. We will 
not take into account the oscillations of the primordial 
spectrum in the rest of the computation and only focus 
on the envelope. The parameter R is closely related to 
the parameter Fb'- high values of R correspond to small 
values of the potential energy of the scalar field at the 
bounce. The link between the phenomenological param- 
eter fco and the physical parameters is less simple since 
the position of the bump depends on both Fb and the 
length of the slow-roll phase. 

The resulting B-mode power spectra are presented as 
a function of the parameter fco in Fig. [9] and [TO] with 
tensor-scalar ratios of 10~^ and 10^^ respectively. At 
large multipole, the spectra are dominated by the con- 
tribution of the lensing as in the usual scenario and are 
thus independant of fco- However, the spectra also reach a 
maximum in two regions: at large angular scales (I < 10) 
for fco < 10^'' and at intermediate scales (20 < i < 200) 
for values of fco around 10~^. The latter case is the most 
interesting from an observational point of view since it 
coincides with the beginning of the observationnal win- 
dow of ground based experiments like BICEP |30j . 

An important question is whether future experiments 
would be able to discriminate between B-modes produced 
by inflation and B-modes arising from a bouncing cosmol- 
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FIG. 10: Angular power spectrum of the B modes as a func- 
tion of fco with a tensor-scalar ratio of 10~^ and R = 500. 
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FIG. 11: Angular power spectrum of the B modes divided by 
the standard inflationary prediction as a function of fco with 
a tensor-scalar ratio of 10~^ and R = 500. 



FIG. 12: Angular power spectrum of the B modes divided by 
the standard inflationary prediction as a function of fco with 
a tensor-scalar ratio of 10~^ and R = 500. 



FIG. 13: B modes for i? = 50 (dashed), R = 500 (long dashed) 
and fco = 10"'^ for a tensor-scalar ratio of 10~^. The solid 
line corresponds to the prediction of inflation with a tensor- 
scalar ratio of 10~^. The 95% confidence upper limits on the 
detection of B-modes of polarization given by BICEP [3D] are 
also pictured. 



ogy in the framework of LQC. A global analysis on all 
multipoles of a similar model was performed in |23j using 
a Fisher analysis method to determine the detectability 
of the bump of the primordial spectrum of gravitational 
waves. Here we will take a complementary point of view, 
we will determine the multipoles where the B-modes aris- 
ing from LQC can be distinguished from those produced 
by a standard inflationary scenario. In Fig. [TT] and [12] 
are presented the B-mode power spectra divided by the 
standard inflationary prediction as a function of the pa- 
rameter fco with tensor-scalar ratios of 10~^ and 10"^ 
respectively. Due to the contribution of the lensing, the 
LQC and inflationary spectra are identical at small angu- 
lar scales. However, for multipoles smaller than 200 the 
LQC spectra show significant departures from the stan- 
dard infiationary prediction for identical tensor-scalar ra- 
tios in the two regions described previously. For R = 500, 
which corresponds to the primordial spectrum of Fig. [6] 
with Fb — 10^^, those departures can reach two orders 
of magnitude for fco ~ 10^^. The corresponding angu- 
lar power spectra are presented in Fig. 13 and 14 In 



the following we will only focus on spectra obtained for 
fco = 10~^ since it corresponds to the region of the pa- 
rameter space where the models can be most easily dis- 
tinguished. 

For high values of R the spectra are very different from 
the inflationary spectra as they possess a peak around 
I = 100 which should allow to discriminate between the 
models. In fact in the case of a tensor-scalar ratio r 
of 10~^, a large region of phase space is even already ex- 
cluded by the BICEP data ^ (see Fig. [T3|. For smaller 
tensor-scalar ratios like r = 10~^ only high values of R 
are excluded by the data. For smaller values of R and 
thus for higher values of Fb the spectra come closer to the 
inflationary prediction in the observable range > 20), 
especially for low values of the tensor-scalar ratio. Nev- 
ertheless, the importance of the backreaction 01] grows 
as Fb gets bigger and the effective model of LQC that we 




FIG. 14: B modes for = 50 (dashed), R = 500 (long dashed) 
and fco = 10~^ for a tensor-scalar ratio of 10~^. The solid 
line corresponds to the prediction of inflation with a tensor- 
scalar ratio of 10~^. The 95% confidence upper limits on the 
detection of B-modes of polarization given by BICEP [30] are 
also pictured. 



are using here is no more valid. Therefore our conclusions 
only concern relatively large values of R. 

The situation changes as the tensor-scalar ratio de- 
creases and the inflationary and LQC spectra might be 
indistinguishable. In Fig. [l5] not only is the gap be- 
tween the inflationary and LQC spectra for r = 10~^ 
much smaller than in the r = 10~^ or r = lO^'^ cases 
but the LQC spectrum for r — 10^^ also yields the same 
B-modes as inflation with r = 10~^ in the observable 
range {I > 20). This arises from the fact that the LQC 
spectra no longer possess a peak at low tensor-scalar ra- 
tios but a plateau around i = 100 which is similar to 
the one present in the inflationary spectrum. The B- 
modes of polarization produced in the frame work of 
LQC should thus be distinguishable in the observable 
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range of multipoles from the inflationary prediction for 
low values of the potential energy at the bounce and for 
moderate tensor-scalar ratios. Therefore, future satellite 
and ground-based experiments should cither be able to 
discriminate between LQC and inflation or to strongly 
constrain the effective model of LQC used in this study. 



FIG. 15: B modes for R = 500 (long dashed) and ko = 10"^ 
for a tensor-scalar ratio of 10~^. The solid and dashed lines 
correspond respectively to the predictions of inflation with 
tensor-scalar ratios of 10~^ and 10~^. The 95% confldence 
upper limits on the detection of B-modes of polarization given 
by BICEP [30] are also pictured. 
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